custEM: Customizable finite-element simulation of complex controlled-source electromagnetic data

We have developed the open-source toolbox custEM (customizable electromagnetic modeling) for the simulation of complex 3D controlled-source electromagnetic (CSEM) problems. It is based on the open-source finite-element library FEniCS, which supports tetrahedral meshes, multiprocessing, higher order polynomials, and anisotropy. We use multiple finite-element approaches to solve the time-harmonic Maxwell equations, which are based on total or secondary electric field and gauged potential formulations. In addition, we develop a secondary magnetic field formulation, showing superior performance if only magnetic fields are required. Using Nédélec basis functions, we robustly incorporate the current density on the edges of the mesh for the total field formulations. The latter enable modeling of CSEM problems taking topography into account. We evaluate semianalytical 1D layered-earth solutions with the pyhed library, supporting arbitrary configurations of dipole or loop sources for secondary field calculations. All system matrices have been modified to be symmetric and solved in parallel with the direct solver MUMPS. Aside from the finite-element kernel, mesh generation, interpolation, and visualization modules have been implemented to simplify and automate the modeling workflow. We prove the capability of custEM, including validation against analytic-solutions, crossvalidation of all implemented approaches, and results for a model with 3D topography with four examples. The object-oriented implementation allows for customizable modifications and additions or to use only submodules designed for special tasks, such as mesh generation or matrix assembly. Therefore, the toolbox is suitable for crossvalidation with other codes and as the basis for developing 3D inversion routines

Evaluation of three approaches for simulating 3D time-domain electromagnetic data

We implemented and compared the implicit Euler time-stepping approach, the inverse Fourier Transform-based approach, and the Rational Arnoldi method for simulating 3D transient electromagnetic data. We utilize the finite-element method with unstructured tetrahedral meshes for the spatial discretization supporting irregular survey geometries and anisotropic material parameters. Both, switch-on and switch-off current waveforms, can be used in combination with direct current solutions of Poisson problems as initial conditions. Moreover, we address important topics such as the incorporation of source currents and opportunities to simulate impulse as well as step response magnetic field data with all approaches for supporting a great variety of applications. Three examples ranging from simple to complex real-world geometries and validations against external codes provide insight into the numerical accuracy, computational performance, and unique characteristics of the three applied methods. We further present an application of logarithmic Fourier transforms to convert transient data into the frequency domain. We made all approaches available in the open-source Python toolbox custEM, which previously supported only frequency-domain electromagnetic data. The object-oriented software implementation is suited for further elaboration on distinct modeling topics and the presented examples can serve for benchmarking other codes

Geophysical analysis of an area affected by subsurface dissolution - case study of an inland salt marsh in northern Thuringia, Germany

The subsurface dissolution of soluble rocks can affect areas over a long period of time and pose a severe hazard. We show the benefits of a combined approach using P-wave and SH-wave reflection seismics, electrical resistivity tomography, transient electromagnetics, and gravimetry for a better understanding of the dissolution process. The study area, "Esperstedter Ried"in northern Thuringia, Germany, located south of the Kyffhäuser hills, is a large inland salt marsh that developed due to dissolution of soluble rocks at approximately 300 m depth. We were able to locate buried dissolution structures and zones, faults and fractures, and potential fluid pathways, aquifers, and aquitards based on seismic and electromagnetic surveys. Further improvement of the model was accomplished by analyzing gravimetry data that indicates dissolution-induced mass movement, as shown by local minima of the Bouguer anomaly for the Esperstedter Ried. Forward modeling of the gravimetry data, in combination with the seismic results, delivered a cross section through the inland salt marsh from north to south. We conclude that tectonic movements during the Tertiary, which led to the uplift of the Kyffhäuser hills and the formation of faults parallel and perpendicular to the low mountain range, were the initial trigger for subsurface dissolution. The faults and the fractured Triassic and lower Tertiary deposits serve as fluid pathways for groundwater to leach the deep Permian Zechstein deposits, since dissolution and erosional processes are more intense near faults. The artesian-confined saltwater rises towards the surface along the faults and fracture networks, and it formed the inland salt marsh over time. In the past, dissolution of the Zechstein formations formed several, now buried, sagging and collapse structures, and, since the entire region is affected by recent sinkhole development, dissolution is still ongoing. From the results of this study, we suggest that the combined geophysical investigation of areas prone to subsurface dissolution can improve the knowledge of control factors, hazardous areas, and thus local dissolution processes

Towards an open-source landscape for 3-D CSEM modelling

Large-scale modelling of 3-D controlled-source electromagnetic (CSEM) surveys used to be feasible only for large companies and research consortia. This has changed over the last few years, and today there exists a selection of different open-source codes available to everyone. Using four different codes in the Python ecosystem, we perform simulations for increasingly complex models in a shallow marine setting. We first verify the computed fields with semi-analytical solutions for a simple layered model. Then we validate the responses of a more complex block model by comparing results obtained from each code. Finally, we compare the responses of a real-world model with results from the industry. On the one hand, these validations show that the open-source codes are able to compute comparable CSEM responses for challenging, large-scale models. On the other hand, they show many general and method-dependent problems that need to be faced for obtaining accurate results. Our comparison includes finite-element and finite-volume codes using structured rectilinear and octree meshes as well as unstructured tetrahedral meshes. Accurate responses can be obtained independently of the chosen method and the chosen mesh type. The runtime and memory requirements vary greatly based on the choice of iterative or direct solvers. However, we have found that much more time was spent on designing the mesh and setting up the simulations than running the actual computation. The challenging task is, irrespective of the chosen code, to appropriately discretize the model. We provide three models, each with their corresponding discretization and responses of four codes, which can be used for validation of new and existing codes. The collaboration of four code maintainers trying to achieve the same task brought in the end all four codes a significant step further. This includes improved meshing and interpolation capabilities, resulting in shorter runtimes for the same accuracy. We hope that these results may be useful for the CSEM community at large and that we can build over time a suite of benchmarks that will help to increase the confidence in existing and new 3-D CSEM codes.Petroleum Engineerin

Capability of low-temperature SQUID for transient electromagnetics under anthropogenic noise conditions

Transient electromagnetics (TEM) is a well-established method for mineral, groundwater, and geothermal exploration. Superconducting quantum interference device (SQUID)-based magnetic-field receivers used for TEM have quantitative advantages and higher sensitivity compared with commonly used induction coils. Special applications are deep soundings with target depths > 1 km and settings with conductive overburden. However, SQUIDs have rarely been applied for TEM measurements in environments with significant anthropogenic noise. We compared a low-temperature SQUID with a commercially available induction coil in an area affected by anthropogenic noise. We acquired four fixed-loop data sets with totally 61 receiver stations close to Bad Frankenhausen, Germany. The high sensitivity of the SQUID enables low noise levels, which lead to longer high-quality transient data compared with the induction coil. The effect of anthropogenic and natural noise sources is more critical for the coil than for the SQUID data. In the vicinity of the transmitter loop, systematic distortion of the coil signals occurs at early times, most probably caused by sferic interferences. We have developed 1D inversion results of both receivers that matched well in general. However, the SQUID-based models were more consistent and showed greater depths of investigation. This led to a superior resolution of deeper layers and even enabled a potential detection of thin conducting targets at up to a 500 m depth. Moreover, we find that the SQUID data inversion revealed multidimensional effects within the conductive overburden. In this regard, we applied forward modeling to analyze systematic differences between inversion results of SQUID and coil data. We determine that low-temperature SQUIDs have the potential to significantly improve the reliability of subsurface models in suburban environments. Nevertheless, we recommend combined application of both types of receivers

A novel semiairborne frequency-domain controlled-source electromagnetic system: Three-dimensional inversion of semiairborne data from the flight experiment over an ancient mining area near Schleiz, Germany

We have developed a novel semiairborne frequency-domain electromagnetic (EM) system and successfully tested it within the DESMEX project. The semiairborne approach relies on the fact that part of the system is positioned on the ground and the rest is airborne. This allows us to take advantage of ground and airborne techniques. In particular, a high-moment transmitter can be installed on the earth's surface, which enables us to inject and induce strong EM fields in the subsurface. Moreover, galvanic coupling is possible, which is an advantage if additional ground stations are deployed. The airborne receivers allow easier, significantly faster, and more uniform spatial coverage of the study area than the ground receivers. In our implementation, transmitters and electric field receivers are installed on the ground. Magnetic field sensors, such as commercially available fluxgate, total field magnetometers, and newly developed induction coils, are installed on a helicopter-towed bird. First, we describe the results of a semiairborne survey performed in a selected area with ancient mining located in the Saxothuringian zone near Schleiz, Germany. A 3D semiairborne inversion model represents several conductive anomalies, which agree well with the outcrop of alum shale formations at the surface. In addition, the shallow parts of the semiairborne model are compared with the result of an independent helicopter-borne survey, which consists of stepwise 1D models

Comparison of novel semi-airborne electromagnetic data with multi-scale geophysical, petrophysical and geological data from Schleiz, Germany

In the framework of the Deep Electromagnetic Sounding for Mineral EXploration (DESMEX) project, we carried out multiple geophysical surveys from regional to local scales in a former mining area in the state of Thuringia, Germany. We prove the applicability of newly developed semi-airborne electromagnetic (EM) systems for mineral exploration by cross-validating inversion results with those of established airborne and ground-based investigation techniques. In addition, supporting petrophysical and geological information to our geophysical measurements allowed the synthesis of all datasets over multiple scales. An initial regional-scale reconnaissance survey was performed with BGR's standard helicopter-borne geophysical system deployed with frequency-domain electromagnetic (HEM), magnetic and radiometric sensors. In addition to geological considerations, the HEM results served as base-line information for the selection of an optimal location for the intermediate-scale semi-airborne EM experiments. The semi-airborne surveys utilized long grounded transmitters and two independent airborne receiver instruments: induction coil magnetometers and SQUID sensors. Due to the limited investigation depth of the HEM method, local-scale electrical resistivity tomography (ERT) and long-offset transient electromagnetic (LOTEM) measurements were carried out on a reference profile, enabling the validation of inversion results at greater depths. The comparison of all inversion results provided a consistent overall resistivity distribution. It further confirmed that both semi-airborne receiver instruments achieve the bandwidth and sensitivity required for the investigation of the resistivity structure down to 1 km depth and therewith the detection of deeply seated earth resources. A 3D geological model, lithological and geophysical borehole logs as well as petrophysical investigations were integrated to interpret of the geophysical results. Distinct highly-conductive anomalies with resistivities of less than 10 Om were identified as alum shales over all scales. Apart from that, the petrophysical investigations exhibited that correlating geophysical and geological information using only one single parameter, such as the electrical resistivity, is hardly possible. Therefore, we developed a first approach based on clustering methods and self-organizing maps (SOMs) that allowed us to assign geological units at the surface to a given combination of geophysical and petrophysical parameters, obtained on different scales. © 2020 The Author

DESMEX: A novel system development for semi-airborne electromagnetic exploration

There is a clear demand to increase detection depths in the context of raw material exploration programs. Semi-airborne electromagnetic (semi-AEM) methods can address these demands by combining the advantages of powerful transmitters deployed on the ground with efficient helicopter-borne mapping of the magnetic field response in the air.The penetration depth can exceed those of classical airborne EM systems,since low frequencies and large transmitter-receiver offsets can be realized in practice. A novel system has been developed that combines high-moment horizontal electric bipole transmitters on the ground with low-noise three-axis induction coilmagnetometers, a three-axis fluxgate magnetometer and a laser gyroinertial measurement unit integrated within a helicopter-towed airborne platform. The attitude data are used to correct the time series for motional noise and subsequently to rotate into an Earth-fixed reference frame. In a second processing step, and as opposed to existing semi-airborne systems, we transform the data into the frequency domain and estimate the complex-valued transfer functions between the received magnetic field components and the synchronously recorded injection current by regression analysis. This approach is similar to the procedure employed in controlled-source EM. For typical source bipole moments of 20-40 kAm and for rectangular current waveforms with a fundamental frequency of about 10 Hz, we can estimate reliable three-component transfer functions in the frequency range from 10-5000 Hz over a measurement area of 4 x 5 km2 for a single source installation. The system has the potential to be used for focused exploration of deep targets

Comparison of novel semi-airborne electromagnetic data with multi-scale geophysical, petrophysical and geological data from Schleiz, Germany

In the framework of the Deep Electromagnetic Sounding for Mineral EXploration (DESMEX) project, we carried out multiple geophysical surveys from regional to local scales in a former mining area in the state of Thuringia, Germany. We prove the applicability of newly developed semi-airborne electromagnetic (EM) systems for mineral exploration by cross-validating inversion results with those of established airborne and ground-based investigation techniques. In addition, supporting petrophysical and geological information to our geophysical measurements allowed the synthesis of all datasets over multiple scales. An initial regional-scale reconnaissance survey was performed with BGR's standard helicopter-borne geophysical system deployed with frequency-domain electromagnetic (HEM), magnetic and radiometric sensors. In addition to geological considerations, the HEM results served as base-line information for the selection of an optimal location for the intermediate-scale semi-airborne EM experiments. The semi-airborne surveys utilized long grounded transmitters and two independent airborne receiver instruments: induction coil magnetometers and SQUID sensors. Due to the limited investigation depth of the HEM method, local-scale electrical resistivity tomography (ERT) and long-offset transient electromagnetic (LOTEM) measurements were carried out on a reference profile, enabling the validation of inversion results at greater depths. The comparison of all inversion results provided a consistent overall resistivity distribution. It further confirmed that both semi-airborne receiver instruments achieve the bandwidth and sensitivity required for the investigation of the resistivity structure down to 1 km depth and therewith the detection of deeply seated earth resources. A 3D geological model, lithological and geophysical borehole logs as well as petrophysical investigations were integrated to interpret of the geophysical results. Distinct highly-conductive anomalies with resistivities of less than 10 Omega rn were identified as alum shales over all scales. Apart from that, the petrophysical investigations exhibited that correlating geophysical and geological information using only one single parameter, such as the electrical resistivity, is hardly possible. Therefore, we developed a first approach based on clustering methods and self-organizing maps (SOMs) that allowed us to assign geological units at the surface to a given combination of geophysical and petrophysical parameters, obtained on different scales. (C) 2020 The Author(s). Published by Elsevier B.V